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Abstract 

We present an algorithm that computes a shortest non-contractible and a shortest non¬ 
separating cycle on an orientable combinatorial surface of bounded genus in 0(n log n) time, 
where n denotes the complexity of the surface. This solves a central open problem in compu¬ 
tational topology, improving upon the current-best 0(n 3 ^ 2 )-time algorithm by Cabello and 
Mohar (ESA 2005). Our algorithm uses universal-cover constructions to find short cycles 
and makes extensive use of existing tools from the field. 


1 Introduction 

Decomposing a surface into topologically simple components is a basic technique in computa¬ 
tional geometry. Many applications in computer graphics, like texture mapping, compression, 
and morphing, rely on decompositions of topological surfaces (see the overview in J21 Chap. 2.4] 
and the references in S]). 

Quite recently, the computational complexity of decomposition problems on combinatorial 
surfaces has gained growing attention. Erickson and Har-Peled 9 investigated the problem 
of cutting a topological surface efficiently into a disk and Colin de Verdiere and Lazarus (Hj 
showed how to decompose a given surface into parts of genus 0. Later works by Erickson and 
Whittlesey m and Cohn de Verdiere and Erickson [1] deal with the computation of short 
generators of the fundamental group and shortest paths and cycles in given homotopy classes. 
A central open question that arose in this context is: 

(P) How to compute a shortest non-trivial cycle on a topological surface efficiently? 

In other words, how to find the smallest cut that reduces the topological complexity of a given 
surface. There actually exist two very natural interpretations of the term “non-trivial,” which 
are both studied in the aforementioned papers: non-contractible and non-separating cycles. 
The former are non-trivial with respect to homotopy, i.e., they are non-zero elements in the 
fundamental group of the surface, and the latter are non-trivial with respect to Z 2 -homology. 

Erickson and Har-Peled [9j presented a first efficient, 0(n 2 log n)-time algorithm for this 
problem on orientable surfaces, where n denotes the size of the surface. It forms a fundamental 
building block, and also the running-time bottleneck, for the path-shortening algorithm in 0J. 
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A natural specialization of the general problem statement is to consider surfaces of fixed or 
bounded genus. This setting excludes the topological complexity from the problem and focuses 
on more essential algorithmic questions. The current-best algorithm for Problem (P) in the 
case of bounded-genus surfaces is due to Cabello and Mohar {2[. They show how to compute 
shortest non-contractible cycles in 0(n 3 / 2 ) time and non-separating cycles in 0(n 3//2 logn) time, 
for orientable and non-orientable surfaces alike. (For the non-separating case, their algorithm 
actually yields subquadratic running time for any g £ 0{n l ^~ e ).) In the full version of their 
paper, they name the task of finding a near-linear-time solution to Problem (P) (possibly for 
surfaces of fixed topological type) as one of the most appealing open questions in the area. 

The present paper solves this problem—on orientable surfaces—for both, the non-contractible 
and the non-separating case. We present an algorithm that computes a shortest non-trivial cy¬ 
cle on a combinatorial surface of bounded genus in O(nlogn) time, where n denotes the total 
description complexity of the surface. 

Our algorithm builds heavily on previous techniques, mainly from ummm- After a first 
introduction of the necessary notions and some formalism, we shall therefore revisit in a little 
greater detail the central algorithmic results from the field. 

2 Topological Background and Formalism 

We begin by fixing our formalism and notations for combinatorial surfaces and thereby also 
briefly recapitulate the most important underlying concepts from topology. This is not meant 
to, and clearly cannot be, an introduction to algebraic topology. So we only go into detail where 
rigorous definitions are indispensable. Beyond that, a fundamental familiarity of the reader with 
the basic concepts of homotopy (and also a little homology) is assumed. We refer to standard 
texts like ,1] or m for a more thorough introduction to the matter. 

A combinatorial surface is a 2-dimensional manifold A4 together with a graph G embedded 
on A4 such that every face (i.e., every component of M \ G) is a topological disk. For surfaces 
with boundary we further demand that every boundary component of A4 be a simple cycle 
in G. All surfaces considered in this paper are implicitly assumed to be connected. 

For a representation of combinatorial surfaces, we follow the formalism of Colin de Verdiere 
and Erickson [1;. We form the dual graph G* of G, which has a vertex /* for every face / 
of (Af,G), and an edge e* for every edge e of G, connecting the two dual vertices ff and /| 
corresponding to the faces f\ and /2 that share the edge e. In case A4 has a boundary, G* also 
contains a vertex e* for every boundary edge e of G, which, in G*, is connected by the edge e* 
to the dual f* of the G-face / that contains e. 

The pair G, G* (together with the correspondences of vertices, edges, and faces of the two 
graphs) contains all topological information about the combinatorial surface (Af,G). As the 
size of the combinatorial surface we define the number of all vertices, edges, and faces of G and 
G* together. In the following, we shall often omit the reference to either the graph G or the 
underlying manifold Ad, simply speaking of “the surface,” implicitly assuming that we have a 
topological space M. represented by a pair G, G* of graphs. 

Usually, the surfaces under consideration in this paper will have no boundary. Only for 
technical reasons we might sometimes have to create boundaries by cutting along cycles. So 
unless stated otherwise, all surfaces are implicitly assumed to come without boundary. 

All the algorithms considered in this paper work for weighted graphs as well as for the 
unit-distance setting. In the former, each edge e of G comes with a non-negative weight, while 
in the latter each edge has unit-weight 1. Edge weights in the dual graph will be defined when 
needed. 
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Paths and cycles on a surface. A path in a manifold Ad is a continuous map a from the 
unit interval [0,1] into Ad. The images of 0 and 1 are considered the start and end point of the 
path, respectively. A loop (with basepoint x) in Ad is a path with a(0) = q(1) = x. A cycle is 
a continuous map 7 from the 1-sphere S l into Ad. So loops and cycles are essentially the same 
thing, only that a loop comes with a distinguished basepoint. One usually identifies a path, 
loop, or cycle with its respective image in Ad. 

In a combinatorial surface, the above objects are embodied by graph-theoretic concepts. A 
path on a combinatorial surface G is a (not necessarily simple) path in the graph G, i.e., a 
sequence of G-adjacent vertices. A cycle is a closed walk. A simple path or cycle contains no 
vertex more than once. The length or weight of a walk is the sum (with multiplicities) of all its 
edge weights. Note that this nomenclature differs from the graph-theoretic convention, where 
a “path” is usually already assumed to be simple and a “non-simple path” would be called a 
“walk.” Here we chose to deviate from that tradition in order to avoid inconsistencies with 
standard topological notions. 

Contractible and separating cycles. Two paths a and /3 are homotopic if there exists 
a continuous map (a homotopy ) transforming a into /?, while keeping start and end points 
fixed. In particular, a and (3 must have the same start and the same end point. Two loops are 
homotopic if they are homotopic as paths and two cycles are homotopic if one can be freely 
transformed into the other, i.e., no points are required to stay fixed. 

The equivalence classes of all loops (with a fixed common basepoint x) form a group— 
the fundamental group 7 Ti(Ad) of Ad—with concatenation of paths as group operation and 
the constant map as identity element. For a connected surface, the fundamental group is 
independent of the choice of the basepoint x. A loop that is homotopy equivalent to the 
constant loop (or cycle) is called trivial (w.r.t. homotopy) or contractible. We call a cycle or 
path tight if it is a shortest one in its homotopy class. 

We distinguish two kinds of non-contractible cycles (or loops) on surfaces: surface sepa¬ 
rating and non-separating cycles. A simple cycle is separating if its removal from the surface 
decomposes the surface into more than one component. Every contractible cycle is separat¬ 
ing. Technically, a separating cycle is trivial with respect to Z 2 - homology. Again, we refer to 
standard text books on topology for an introduction to homotopy and homology. 

The genus of a surface. By now, our results seem to work only for orientable manifolds, so 
we assume from now on that all surfaces under consideration are orientable. The topological 
type (up to homeomorphism) of an orientable surface Ad is uniquely determined by its genus 
g > 0 and (in case Ad has a boundary) by the number b of boundary components. Intuitively, 
g counts the number of “holes” in the surface. For example, the torus has genus 1. 

It is well-known that the fundamental group of a genus -5 surface without boundary is gener¬ 
ated by 2 g loops and moreover, one can choose the loops in such a way that each of them forms a 
simple non-separating cycle and any two of them intersect only in the common basepoint. If we 
now cut along all these loops, the result will be a topological disk. Erickson and Whittlesey m 
call such a collection of 2 g almost-disjoint loops that generate 7 Ti(Ad) and whose removal leaves 
a disk, a system of loops. Figure ^ shows a system of loops on a genus-2 surface. 

Technically, the existence of a system of loops is guaranteed only for topological spaces; on a 
combinatorial surface, there might simply not be enough space to route 2 g loops without large 
overlaps. However, it turns out that these issues can be overcome. We shall address this point 
in more detail later. 
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Figure 1: A system of loops on a 2-hole torus. 


The universal cover. The universal cover of a connected topological space JA is the (up 
to homeomorphism) uniquely determined simply-connected space JA such that there exists a 
map / from JA onto JA such that / is locally a homeomorphism. The universal cover of any 
orientable surface of positive genus (and without boundary) is the infinite plane. 

Fix a system of loops l\,...,l 2 g for JA. The universal cover JA is tiled with infinitely 
many copies of JA — (the surface JA cut along the loops If), each such copy D being 

called a. fundamental domain. Two adjacent fundamental domains touch along a loop Every 
fundamental domain has each loop twice on its boundary, seeing it once from each side. If we 
orient each loop (arbitrarily) then its “left” and “right side” are well-defined and we denote 
their two occurrences on the boundary of a fundamental domain D by if and £~, depending 
on whether D lies to the left or right of the respective copy of i. See Figure [3 



Figure 2: A fundamental domain of a 2-hole torus. 

A path a in JA can be lifted to a path in the covering space by choosing any inverse image 
of the starting point of a and then “unrolling” a into JA, making a transition between two 
adjacent fundamental domains whenever the path in JA crosses a loop from the respective 
system of loops. This gives a procedure to decide whether two paths a and /3 in JA (with the 
same starting and end point) are homotopy-equivalent: After choosing the same starting point 
for a lift of a and /? into JA, we see that a and f3 are equivalent if and only if their lifts also 
end on the same point in JA. In particular, a loop in JA is contractible if and only if any, and 
thus all of its lifts are also loops. 

3 Previous Algorithms 

During the last few years, some powerful tools for the computation with paths and cycles 
on combinatorial surfaces have been developed. Several of those techniques will be used in our 
algorithm. Let us first briefly sketch the relevant results to introduce the reader to the methods. 
Later, when we present our new algorithm, we will address more technical details. 
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Shortest non-trivial loops with given basepoint. In [S]. Erickson and Har-Peled de¬ 
scribed a fundamental method for finding a shortest non-contractible cycle through a given 
vertex x on a combinatorial surface. Essentially, they just run Dijkstra’s shortest-path algo¬ 
rithm, with starting point x. Whenever the wave front of equidistant vertices touches itself as 
it sweeps across the surface, the algorithm checks whether it does so in a trivial way, enclosing 
a topological disk on either side, or whether it found two homotopically different paths to that 
contact point, in which case a non-contractible loop has been found. 

Altogether, a shortest non-contractible cycle containing a given vertex can be found in 
0(n logn) time (where n denotes the complexity of the surface). With slight modifications, a 
shortest non-separating cycle can be found in the same time. The authors remark that this 
immediately gives an 0(n 2 log n) algorithm for finding a shortest non-trivial cycle by simply 
running the above algorithm for once for each vertex of the surface. 

Constructing a whole system of loops. The above algorithm for shortest-loop finding can 
be used to compute a whole system of loops. After computing a shortest non-separating loop 
l\ through x, one can cut along this loop (implicitly duplicating every vertex on i\) and then 
compute another “loop” from one copy of x to the other in the new surface Ad — i\. Intuitively, 
the cut along t\ forces the next loop, £ 2 , to lie in a different homotopy class than l\. We may 
repeat this process of finding loops and cutting until 2 g loops are found, which can easily be 
seen to generate the fundamental group of the surface. The resulting loops will in general not 
be disjoint but this technical detail is not too hard to resolve. 

Erickson and Whittlesey m prove the remarkable fact that the greedy construction above 
always produces a shortest system of loops (i.e., a system of loops that minimizes the total length 
of all loops) for the given basepoint. They also present an alternative approach for computing 
a shortest system of loops through a given vertex in linear time for a surface of bounded genus. 
That construction uses Eppstein’s tree-cotree decomposition, which we recapitulate next. 

The tree-cotree decomposition. In [J], Eppstein presented a dynamic data structure for 
maintaining the embedding of a graph on a surface, which turned out to be a very powerful tool 
even in the static setting. For a given surface At, two trees are computed: a rooted spanning 
tree T (with root x) for the graph G and a spanning tree T* for the dual graph G*. The two 
trees must not intersect, in the sense that if an edge e occurs in T then its dual edge e* must 
not occur in T*. and vice versa. 

Let X denote the set of all edges of G that do not occur in T and whose duals are not part 
of T*. By Euler’s formula, |Aj = 2 g. The triple ( T,T*,X ) is called a tree-cotree decomposition 
for G. Let P(u) denote the unique path in T from a vertex u to the root x of T. Each edge 
e = {it, v} € X induces a loop £ e = P(u) U e U P(v) in T. Eppstein proves that these 2 g loops 
form a homotopy basis and shows that for weighted graphs with all edge weights distinct, a 
minimum-spanning tree T in G and a maximum spanning tree T* in the dual G* are disjoint 
(in the above sense that at most one of e and e* occurs in either tree) and can thus serve as 
tree and cotree, respectively, to form such a decomposition. 

Erickson and Whittlesey m use the tree-cotree decomposition to obtain a shortest system 
of loops for a given basepoint x very efficiently. First they compute a shortest-paths tree T, as 
above. Then they assign new weights w(e) to the remaining edges e £ E(G) \ E(T ): for such 
an edge e let u(e) be the total length of the induced loop £ e . They then compute a maximum 
spanning tree T* in the dual (G \ T)*, using as the weight of a dual edge e* the weight uj(e) of 
e, too. Again, an edge set X of size 2 g will be left and it can be shown that now {£ e \ e € X} 
will actually be a shortest system of loops. (Again, the issue of overlapping loops is deferred.) 
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A great advantage of the tree-cotree decomposition over a repeated application of Erickson 
and Har-Peled’s non-trivial-loop algorithm is that it runs in linear time on surfaces of bounded 
genus. More precisely, Erickson and Whittlesey observe that using an algorithm of Henzinger et 
al. EH, shortest-paths trees can be computed in linear time on a surface with genus g € (^(n 1 e ). 
Note that “computing” here means providing the tree T together with the set X that induces 
the 2 g loops. Simply outputting all resulting loops explicitly might require already £l(gn) time. 

Routing through finite portions of the universal cover. A standard technique for find¬ 
ing non-contractible cycles is to compute shortest paths in a finite portion of the universal cover. 
Assume we have a collection of cycles 7 \ whose removal decomposes our surface into a collection 
of topological disks. A system of loops, for example, has this property (producing exactly one 
disk). In the universal cover At, each cycle C{ appears as a collection of infinite paths. If we 
now have some bound on the number of times a potential path or cycle we are looking for can 
cross the cycles Ci, we know that a lift of the desired path or cycle can be found in a finite part 
of the universal cover. 

Cabello and Mohar f2J use such an approach to find shortest non-contractible cycles and 
Colin de Verdiere and Erickson Q compute shortest paths and cycles within a given homotopy 
class by similar methods. While Cabello and Mohar work with systems of loops, Colin de 
Verdiere and Erickson perform a so-called tight octagonal decomposition of the given surface such 
that the universal cover becomes a union of octagonal regions that are separated by shortest 
paths. 

4 Our Algorithm for Shortest Non-Trivial Cycles 

Our algorithm starts by computing a shortest system of loops as demonstrated in m- Precisely, 
we first compute a shortest-paths tree T from some arbitrarily chosen basepoint x and then find 
a maximum-spanning tree T* for the remaining dual graph with edge weights c o as discussed 
in the previous section. This takes 0(n) time if the genus g of the surface is 0(n 1 ~ e ) m 
Thru. 3.9]. 

Let us address the case of overlapping loops in a bit more detail. Cabello and Mohar j2|, 
who also use the tree-cotree decomposition to find short cycles, describe a simple transformation 
of the underlying graph G that turns all loops into simple cycles (still with basepoint x ) that 
pairwise only intersect in x. 

By construction, the two paths P{u) and P(v) that make up the loop t e for an edge e = { u , v} 
from the set X , defined as above, intersect only on a path a starting from the basepoint x until 
they follow different branches of the tree T. Cabello and Mohar duplicate all vertices along 
a (except x) introducing zero-weight edges between corresponding vertices so that afterwards 
P(u) and P{v) only touch at x, and i e thus becomes a simple loop. Similarly, they duplicate 
vertices wherever different loops £ e and £ e / for e, e' E X intersect. This whole process introduces 
no more than 0(gn ) new edges and vertices and also requires 0(gn ) time 0 Lem. 7]. 

One easily sees that the “split” loops l e (e £ X) now form a system of loops for the modified 
graph, which clearly maintained the topology of the original one. In fact, in our situation, this 
will still be a shortest system of loops because we have only introduced zero-weight edges into G. 
Formally this can be seen by running the whole construction of the tree-cotree decomposition 
from scratch for the modified graph. Taking care that none of the new edges is used for the 
tree T (which can easily be guaranteed) we reproduce exactly the same system of loops that we 
obtained through the splitting process—which must now be a shortest system of loops by m 
Sec. 3.5]. 
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Controlling crossings. A central ingredient for the shortest non-contractible-cycle algorithm 
in J2j is Lemma 7 from that paper, which guarantees the existence of a shortest non-contractible 
cycle that intersects the constructed system of loops not too often. We will employ that result 
for our algorithm, too. However, our formalism will deviate slightly from the one in j2] because 
we have to prepare for later needs and we also have to extend that lemma to non-separating 
cycles (which was implicitly mentioned in [2j already). Moreover, the respective lemma was 
stated with a very sketchy proof only. Therefore, we will present it here again, with full proof. 

Intuitively it is clear what a “crossing” between two cycles on a surface should be. Yet, 
formally, we have to make precise whether and how two cycles are considered to “cross” if they 
share several consecutive edges. To this end, we think of every loop £i (and thus every of its 
edges) to have two “sides” and a path or cycle that uses such an edge e from one of these 
loops is considered to he by an infinitesimal margin to the left or right of this edge e. In the 
following, all cycles are implicitly assumed to come with fixed side informations for every single 
edge of them that lies on one of the loops Figure 01 shows a figurative drawing of this extra 
information. 



Figure 3: Every cycle contains side information for edges on the loops £{. 

Now, a crossing between a cycle 7 and a loop ti is a pair of consecutive edges e, / in 7 
such that the vertex of 7 at which e and / touch lies on ti and such that locally, e and / lie 
on different sides of £i, taking the infinitesimal offsets into account, of course. The path in 
FigureOl for example, crosses ii exactly twice, £j just once (at x ), and £h three times. 

There is a problem with cycles through the basepoint x. Since the loops £ t only intersect at 
this point, this is the only vertex at which some cycle could intersect more than one loop at the 
same time. Such a situation would bring about some inconvenient technicalities, so we modify 
our graph G a bit now in order to be able to ignore this special case in the future. 

We subdivide every edge e incident with x into two edges and then connect the subdivision 
points in a cycle, in the order they occur on the surface. The weights of all the edges that are 
now incident with x are set to some large number L greater than all other edge weights in G 
together, while the edges on the new cycle receive zero weight; the outer edges maintain the 
weight of the respective original edges. See Figure El 

This modification obviously does not change the topology of the surface and it also does not 
change the distances between vertices—except those from and to x, of course. The advantage of 
the modified graph is that while the loops £ 1 , ■■■ ,£ 2 g still all go through the origin, any shortest 
(non-trivial) cycle on the surface will now avoid the origin. This is easy to see. If some cycle 
went through x in the original graph, we can reroute it on the zero-length cycle around x now, 
to obtain a homotopically equivalent cycle of the same length as the original one. Conversely, 
any cycle in the modified graph maps to an equivalent one of same length in the old graph by 
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Figure 4: Blowing the basepoint up to a disk. The labels show edge weights. 


shifting all vertices on the new cycle to x and deleting repetitions. Our system of loops will still 
originate at x and it will still be a shortest system of loops. 

The following lemma is a direct adaptation of Lemma 8 from f2_. 

Lemma 1. Let li,...,l 2 g be a system of loops for a surface M as computed above. Then 
there exists a shortest non-contractible and a shortest non-separating cycle on M., each of them 
crossing any loop i{ at most twice. 

Proof. We treat the non-contractible and the non-separating case simultaneously, speaking of 
“non-trivial” cycles to mean both variants and addressing either type explicitly when necessary. 
Note that we only need to consider simple cycles because any self-intersecting non-trivial cycle 
contains a proper subpath that forms a non-trivial cycle. 

Amongst all shortest non-trivial cycles consider one, 7 , that minimizes the total number of 
crossings with all li. If 7 intersects each of the loops at most twice then we are done. Otherwise, 
consider a loop lj that has at least three crossings with 7 . Two of these crossings must lie on 
the same shortest path P{a) in the tree T between x and one vertex a of the edge that defined 
the loop lj. Denote the P(a)-vertices at the two crossings of 7 with lj by u and v, denote the 
subpath of P(a) between u and v by a, and call the two “halves” of 7 (separated by u and v) 
71 and 72 . See Figure El 



Figure 5: The situation from the proof of Lemma ^ 

We have |a| < I 71 1, ( 72 1 because a as a subpath of P(a) is a shortest path. 

We claim that not both of the composite cycles 071 and 072 can be trivial. Assume for 
contradiction that they are. If they are both contractible then their composition 7 iaaT 1 7 ^~ 1 , 
which yields the original cycle 7 , must be contractible, too—a contradiction. In case the two 
cycles are separating, their sum (in Z 2 -homology) also gives 7 , and thus a contradiction as well. 

So we know that one of the two cycles 0:71 and 072 is non-trivial. Assume w.l.o.g. that the 
former is. This cycle is no longer than 7 and we can choose the side of the lj edges it uses in 
such a way that it has fewer crossings with lj than 7 . If 71 connects to a at a and v from the 
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same side of P(a) then we place our new cycle completely to that side of lj, thereby removing 
the two crossings at u and v and introducing no new ones. If 71 connects to a from different 
sides, we route the new cycle along an arbitrary side of lj, which results in just one crossing, 
instead of two. In both situations, we do not create any new crossings with the other loops 
£i because the path a is completely disjoint from them. (At this point we make use of our 
modification of G around x, which guarantees us that 7 does not touch x .) □ 

Lemma [2 tells us that when searching for a shortest non-trivial cycle, we need only consider 
such cycles that cross the system of loops a bounded number of times. Set in the universal cover, 
this means that we need only consider paths that pass through at most 4g + 1 fundamental 
domains there. Naively applied, this approach would be quite inefficient, though. 

Note that the above lemma used only the fact that T was a shortest-paths tree, not that 
the cotree T* was chosen in any particular way. In the following we shall exploit the additional 
fact that the cotree T* was constructed in such a way that the resulting system of loops is a 
shortest one. 

Shortest cycles around cylinders. In Hj, Colin de Verdiere and Erickson introduced a very 
useful lemma for tightening cycles. Using a planar s-t cut algorithm by Frederickson EH, which 
improved upon techniques by Reif El, the authors explain how to compute a shortest cycle on 
a cylinder S 1 x [ 0 , 1 ] in almost linear time. 

Lemma 2 (Colin de Verdiere & Erickson [4, Lem. 3.5(d)]). For a cylinder S 1 x [0,1] of 

(combinatorial) size n, we can compute a tight cycle homotopic to the two boundaries S 1 x {0} 
and S 1 x { 1 } in 0{n log n) time. 

Similar to [1], we will use this cylinder lemma as a tool to find short paths in the universal 
cover. However, we need a little more preparation before we can apply it to our situation. 

Lifting to unique fundamental domains. Let us make the intuitive notion of path finding 
in the universal cover more precise. If we cut the graph G along our system of loops, we 
obtain a plane graph G that has every loop vertex and edge twice on its boundary, except for 
x which occurs exactly 4 g times (see Figure |2] again). We envision the universal cover M of 
J\4 as represented by an infinite plane graph G, which is made up of infinitely many copies of 
fundamental domains G, glued together along copies of the loops li. 

Consider a lift a of some path or cycle a in G that passes through several fundamental 
domains in G. Those edges of a that do not lie on one of the loops ii lie in a uniquely 
determined fundamental domain. For edges along the loops there are in principle always two 
possible fundamental domains they can be assigned to. However, the side information of a 
resolves this ambiguity; in the obvious way: an d-edge on the left of a loop belongs to the 
fundamental domain to the left and an edge to the right belongs to the right-hand-side domain. 
This way, every lift a travels through a uniquely determined sequence of fundamental domains 
in the universal cover G. If a also avoids the basepoint x (as is guaranteed with our graph) 
then any two fundamental domains in such a sequence are adjacent, i.e., they touch along the 
lift of some loop lj. 

Definition 1. We call a sequence V = (Dq, .. . ,D^) of fundamental domains Di in G a meta 
path if every pair Dj_i, Di is adjacent (i.e., touches along a lift of some loop lj). Let a be any 
path or cycle in G. A meta path T> is a meta lift of a if some lift of a travels through the 
fundamental domains in T> in the respective order, i.e., it starts in Dq, then enters D 1 , then D 2 , 
and so on, until it reaches its end point in Dk- 
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Forming cylinders. Algorithmically, we want to use meta paths to find cycles with the help 
of Lemma |2I Therefore, we have to turn meta paths into cylinders. 

The key idea is quite simple: identify the domains Dq and Dj, of some meta path—in the 
right way. Figure El shows a meta path through the universal cover of a torus, for which this 
merging can be done without any problems. Note that in both drawings, D\ borders on the top 
loop of Dq and Dq on the right loop of Da- As we shall see, cycles around such a cylinder will 
produce non-trivial cycles on the underlying torus. 


D a 

d 3 



d 2 

Da 



Do 



Figure 6: Merging the ends of a meta path. 

Several details have to be taken into account for this gluing step, however. First of all, 
we have to deal with the possibility that some of the fundamental domains along V could 
be identical. In such a situation we always create a new copy of the fundamental domain so 
that Dq,D\, ... ,Dk become all distinct. Technically, we do the following. For a given meta 
path D = (Dq, ..., Dk) we take pairwise disjoint copies D ' 0 ,..., D' k of the fundamental domains 
Dq, ..., Dk and glue every pair D(_ 1 , D( along the boundary loops at which the corresponding 
fundamental domains Di_i, Di touch, i.e., they get connected just like in the universal cover. 
For the constructed space, which is obviously a topological disk, we eventually identify the 
domains D' 0 and D' k . Denote the resulting space by T>°. 

Is this space V °always a cylinder? Unfortunately not. If some path a crosses the same 
loop ij twice, first from one side and then from the other, without crossing any other loop 
in between, then the respective steps in a meta path of a will be a sequence ZU_i, Di, Dj+i 
with = D{ + 1 . In other words, a lift of a would enter some fundamental domain and then 

leave it through the same boundary of D j. The resulting space V° would thus contain the 
non-manifold configuration of Figure [3 



Figure 7: A non-manifold situation in the space T>°. 

However, we can show that there will always be a shortest non-trivial cycle that does not 
“curl” into a fundamental domain as the one in Figure 0 and will thus yield a space T>° that is 
a manifold. 
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Definition 2. We call a non-contractible cycle 7 in G curl-free if for every meta lift T> = 
(Do,..., Dk) of 7 we have D*_ 1 / Dj + 1 for 1 < i < k. 

Lemma 3. Amongst all tight cycles in the homotopy class of some given non-contractible cycle, 
the ones that cross the system of loops the least number of times are curl-free. 

Proof. Consider a tight non-contractible cycle 7 that is a shortest cycle in its homotopy class and 
which amongst all such tight cycles intersects the system of loops a minimum number of times. 
Assume for contradiction that for some lift 7 of 7 the corresponding meta lift T> = (Do,..., Dk) 
has an index i with D,_ 1 = Dj + 1 . Let ip denote the subpath of 7 between the two crossing 
points u and v where 7 enters and leaves, respectively, the fundamental domain D*. So <f> “curls” 
into Dj, entering and leaving that fundamental domain through the same boundary loop, tf or 

ij, sa y- 

Replace ip in 7 by the subpath of lj between u and v, which cannot be longer than ip because 
lj is from a shortest system of loops and the two paths are homotopic uni Lem. 3.8]. Routing 
the new path segment on the D,_i-side of lj, we ensure that the resulting path j 1 now does 
not enter Dj at all and thus has two fewer crossings with lj than 7 . A contradiction to the 
minimality assumptions for 7 □ 

Lemma |21 can be seen as an extension of Lemma |U Now we know that some shortest non¬ 
trivial cycle in G will not only cross each loop li at most twice; its lift in G will also always 
leave a fundamental domain through a different loop than the one through which it enters. 

Note that this observation holds for non -separating cycles, too, not only for non-contractible 
ones, although Lemma [21 is only about the homotopy case. This is because we first apply 
Lemma|T]to learn that some shortest non-trivial cycle 7 intersects every loop at most twice and 
then conclude that in the homotopy class of 7 a shortest cycle with the minimum number of 
crossings must be curl-free, irrespective of whether 7 was separating or not. 

Before we give a formal proof that for some shortest non-trivial cycle the space T>° will 
indeed be a cylinder, let us slightly modify that space to make the proof a little easier. Recall 
that by our modification around the basepoint x , no shortest non-trivial cycle will pass through 
x anymore. So we always delete all copies of x, together with the whole disk of incident faces, 
from the space T>°. This modification has no effect on the eventual algorithm but it makes the 
proof of the following lemma a bit easier. 

Lemma 4. There exists a shortest non-contractible (non-separating) cycle 7 with some meta 
lift V = (Do,..., Dk) with 2 < k < 4g such that the space T>° forms a cylinder. 

Proof. Lemma Q guarantees the existence of a shortest non-trivial cycle in G that crosses each 
loop li at most twice and Lemma |21 tells us that if amongst all such cycles we take one, 7 , that 
crosses the system of loops the least number of times then it will also be curl-free. 

Consider a meta path D = (Do,... , D&) of 7 By the crossing bound above we know that 
k < 4g and a lift of any non-trivial cycle must pass through at least two different fundamental 
domains, hence k > 2. The curl-freeness of 7 guarantees that T>° is a manifold. There are only 
three possible topological types that could in principle result from the gluing process: a cylinder, 
a Mobius strip, or a disk. The orientability of A4 excludes a Mobius strip and the deletion of 
all copies of the basepoint implies that the resulting space cannot be simply connected. Hence, 
it must be a cylinder as claimed. □ 

The next lemma shows that the cycles we want to compute in the meta cylinders T>° corre¬ 
spond to short cycles in G as intended. 
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Lemma 5. Let a be a curl-free non-contractible cycle in G and let V be a meta lift of a. Then 
the lift of a into D forms a cycle (3 in V° homotopic to the two boundaries of V° and of the 
same length as a. Conversely, every non-contractible cycle (3 in V° homotopic to the boundaries 
of the cylinder corresponds to a cycle of the same length in G and homotopic to a. 

Proof (sketch). We can cut any cycle (3 in D° at an arbitrary point in Dq to obtain a path 
/ 3' from a point bo in Do to a point b/. in D through the universal cover. Since bo and b *. 
correspond to the same point in A4, such a path /3 r corresponds to a cycle in M. Obviously, 
homotopic cycles in D° correspond to homotopic cycles in At. □ 

The algorithmic realization of the above results is now straight-forward. List all different 
meta paths of length at most 4 g in the universal cover G of G. Form the space T>° from each 
of those meta paths and throw away all those spaces that are not cylinders. Use Lemma [2] 
to compute a shortest non-contractible cycle around each cylinder T>°. This way we obtain a 
collection of cycles 7 in G. By Lemma 03 we know that one of these cycles must be a shortest 
non-trivial cycle of G. 

Since a shortest non-trivial cycle must be simple, we throw away all cycles with self¬ 
intersections. For each of the remaining cycles, we can determine in linear time whether it 
is contractible or separating HI Lem. 5.1 ]. 1 This way we are bound to find a shortest non- 
contractible and a shortest non-separating cycle. 

Theorem 1. For an orientable combinatorial surface of bounded genus, a shortest non-con¬ 
tractible cycle and a shortest non-separating cycle can be found in 0(n log n) time. 

Proof. We have already argued about the correctness of the above algorithm. The running-time 
bounds are easily checked. Our algorithm considers 0((4g) 49 ) meta paths, which is a constant 
if the genus is bounded. A single cylinder computation costs 0{gn log gn) time by Lemma |3 
so the total running-time, for bounded g, is also 0(nlogn). The tests for contractible and 
separating cycles takes linear time per cycle. □ 

Surfaces with Boundary. We can use an idea from (Hj to extend Theorem El to surfaces 
with boundary. In a preprocessing step, we glue a torus into each boundary component of 
M. Precisely, for a boundary cycle C with r edges, we build a genus-1 surface with a length-r 
boundary and attach it along C. All edges inside the new surface parts receive a common weight 
of some number L larger than the sum of all edge weights in the original graph G. The whole 
process, for all boundary cycles together, will increase the combinatorial size of A4 by a linear 
factor, only. It is easy to see that a shortest non-trivial cycle on the modified surface will also 
be a shortest such cycle on the original surface and vice versa. 

An application: tightening paths and cycles. In 0, Colin de Verdiere and Erickson 
present algorithms to find shortest paths and cycles in given homotopy classes of an orientable 
surface. After precomputing a special decomposition of the surface into octagonal slabs in 
0 (n 2 log n) time, they can find a shortest path homotopic to a given path of complexity (combi¬ 
natorial length) k in O(gnk) time and a shortest cycle homotopic to a given cycle of complexity 
k in 0(gnklog(gnk)) time. 

We observe here that Theorem El can be used to speed up the precomputation of their 
algorithm from 0{n 2 logn) to 0(n log n), so that the running time of this step is now dominated 
by the later path computation and can thus be ignored. 

x One could actually just read the relevant information directly off the meta paths but we prefer to avoid the 
necessary group theoretic argument here. 
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5 Conclusion 


We have shown how to compute shortest non-contractible and shortest non-separating cycles 
on orientable surfaces of bounded genus in 0{n log to) time. A question that naturally arises 
from this result is whether it can be translated to the case of non-orientable surfaces. Another, 
possibly more challenging problem would be to bring down the exponential dependence of the 
running-time on the genus of the surface. Maybe a clever inspection of the universal cover might 
yield a polynomial dependence and thus generalize the result to useful statements for surfaces 
of unbounded genus. 
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